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An  air-sea  coupled  model  developed  in  this  paper  is  intended  to  depict  the 
generation  of  unstable  coastal  water  waves  near  a  fast  shelf-ice  edge  with  small 
curvature.  Without  any  external  forcing  of  the  system,  the  unstable  modes  are 
excited  by  the  air-sea  feedback  mechanism.  Thermally  forced  surface  winds, 
generated  by  a  surface  temperature  gradient  that  is  related  to  the  vertical 
displacement  of  the  pycnocline,  in  turn  produce  a  surface  water  current 
that  further  changes  the  vertical  displacement  of  the  pycnocline.  The  model 
consists  of  two  elements:  thermally  forced  surface  air  flow  and  mechanically 
driven  coastal  water  waves  (Clarke).  The  two  components  are  linked  through  the 
surface  temperature  gradient  and  surface  wind  stress.  The  coupled  air-sea  model 
is  solved  as  an  eigenvalue  problem.  It  is  found  that  the  iceward  convex  ice-edge, 
relatively  thin  ice,  and  small  upper-layer  water  depth  favor  the  generation  of  very 
unstable  modes.  The  influence  of  the  coupling  on  the  wave  velocity  is 
also  discussed. 


|.  INTRODUCTION 

' '  I 

j  Clarke  (1978)  discussed  in  detail  the  wind-driven  quasi-geostrophic  water 
movement  near  a  fast  shelf-ice,  i.e.,  glacial  ice  which  is  attached  to  a  land  mass. 
Water  motion  near  the  shelf-ice  edge  is  generated  not  only  by  a  stress  curl  mech¬ 
anism,  but  also  by  a  blocking  action  as  in  the  coastal  case.  The  quasi-geostrophic 
motion  can  be  described  by  wind-forced  trapped  long  waves  which  propagate 
along  the  ice-edge  with  the  ice  on  their  left  (right)  in  the  Southern  (Northern) 
Hemisphere  (Clarke,  1978).  In  Clarke's  model  the  wind  stress  is  taken  as  a  known 
function  (external  forcing).  However,  in  this  paper  the  surface  winds  are  deter- 
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FIGURE  1.  Diagrammatic  representation 
of  two-layer  shelf-ice  model  (Clarke,  1978). 


mined  by  the  surface  temperature  gradient  which  is  related  to  the  vertical  dis¬ 
placement  of  the  pycnocline,  h.  The  air-sea  feedback  model  depicted  in  the 
subsequent  sections  is  intended  to  change  the  Clarke  equation  into  a  homoge¬ 
neous  equation  for  h,  and  then  treat  it  as  an  eigenvalue  problem. 


II.  AIR-SEA  COUPLED  SYSTEM 

Clarke  (1978)  developed  a  two  layer  shelf-ice  model  to  discuss  wind-driven 
quasi-geostrophic  water  movement  near  fast  shelf-ice  edges  with  small  curvature 
(Fig.  1).  The  small  curvature  means  that  the  radius  of  curvature  of  ice-edge  at  any 
point  is  much  greater  than  the  baroclinic  radius  of  deformation.  Such  curvature 
effects  can  be  examined  by  using  normal  and  tangential  coordinates  (n,  s).  Here, 
n  is  defined  as  the  distance  seaward  from  the  ice-edge,  and  s  is  the  distance  along 
the  ice-edge  from  the  origin.  The  curvature  n(n )  is  defined  to  be  positive  (neg¬ 
ative)  when  the  ice  edge  is  convex  seaward  (iceward).  For  small  curvature  of 
ice  edges  the  upward  displacement  of  the  pycnocline  h*  is 


h*(n,  s,  t)  ■  h(s,t)  e 


-  n/»(X  +  k/2) 


♦  0(<2), 


(1) 


where  h(s,  t)  satisfies 


3h  Bh  y(p  -1)  +  3t 

T“  "  *7”  "  - u — i  [  fe  *x  +  e  •  —  ], 
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with 
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where  Hl(  H2  are  the  upper-  and  lower-layer  depths;  a  is  the  radius  of  defor¬ 
mation  of  the  baroplinic  mode;  H2  is  the  lower-layer  depth  under  the  ice;  H  is 
the  total  water  depth;  AH2  =  H,  —  H2 ,  and  AH  =  H  —  H2 ;  k  is  the  ice-edge  cur¬ 
vature  non-dimensionalized  with  respect  to  the  reciprocal  of  the  baroclinic 
radious  of  deformation;  r  is  the  wind  stress;  en  and  e”s  are  the  unit  vectors  in  the 
direction  of  increasing  n  and  s,  respectively;  g'  is  reduced  gravity  g(p2  -p,)/p2. 
Here  g  is  acceleration  due  to  gravity,  and  p,  and  p2  are  the  upper-  and  lower-layer 
densities. 

The  sea  surface  temperature  (SST)  anomaly  Tj  caused  by  pycnocline  fluc¬ 
tuation  is  assumed  to  be  proportional  to  h*  in  the  following  form 


T 

s 


t 


AT  AT  +  k/2) 

—  h*  ■  —  he 


(5) 


which  means  that  the  coastal  upwelling  {h  >  0)  near  a  fast  shelf-ice  will  give  rise 
to  the  warrfi  anomaly.  Here  AT  =  T2  —  T,  >0,  is  a  temperature  difference  be¬ 
tween  the  lower-  and  upper-layers.  Such  a  differential  SST  anomaly  generates  a 
thermally  forced  wind.  Chu  (1986abc)  developed  a  coupled  air-ice-ocean  model 
for  the  marginal  ice  zone  (MIZ)  to  show  that  thermally  forced  surface  winds  blow 


Northern  H«ml*pher»  Southern  Hemisphere 

FIGURE  2.  Relationship  between  thermally 
forced  surface  wind  and  the  surface  tempera¬ 
ture  gradient. 
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across  the  isotherms  from  cold  to  warm  regions  with  some  deflection  due  to  the 
Earth  rotation.  The  deflection  angle  a  (Fig.  2)  depends  on  the  atmospheric 
Richardson  number  (Chu,  1986c),  and  the  wind  speed  is  proportional  to  ivT's i, 
which  is  computed  by 

ivVi  2l(  j  +i)  a-— > .  (6) 

3n  a  2 

Therefore,  the  thermally  forced  surface  winds  are  estimated  by 


.  -  n/a(l  +  k/2) 

a  h  e 
s 


The  parameters  as  and  an  are  defined  by 


anh  e 


•  n/a(l  +  k/2) 


(7) 


a 
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(1  +  */2) 

-  sin  a  , 

a 
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AT  (1  +  k/2) 

a  s  b  —  — - cos  a, 

n  H1  a 


(9) 


where  b  is  a  proportionality  constant  which  represents  the  wind  speed  thermally 
forced  by  a  unity  of  surface  temperature  gradient.  According  to  Defant  (1950) 

1  C/100  km  surface  temperature  gradient  can  generate  about  2  m/s  surface  wind 
Hence,  we  take 


2  m/s 

Pc/lOO  km 


Surface  wind  stress  r  is  then  obtained  by 


T 

S 


P  C  v(a), 

a  a  s 


T 

n 


p  C  v(a) 
a  a  n  1 


(10) 


where  Ca  is  the  dimensional  (m/s)  air  drag  coefficient.  Substituting  (10)  into  (2) 
and  using  (7)  lead  to  the  .following  equation 
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where 


u)  = 


P  c 

Ka  a 
P1H1 


The  amplitude  of  the  pycnocline  displacement  h  is  then  governed  by  a  free, 
first-order  wave  equation  with  two  air-sea  coupling  coefficients  as  and  an. 


I.  WAVE  SOLUTION 

Let  the  solution  of  the  basic  equation  (11)  have  the  following  wave  form: 


hCs,  t)  -  A  e1(ot  ‘  “>• 

Substituting  (12)  into  (11)  gives  the  dispersion  relation 


YPg  m 


fy(y  -  l)w  a 


Wave  speed  is  computed  by 


Re(o) 


m 


Growth  rate  is  given  by 


YPg 


-  Im(o)  ■ 


yg'  +  y(p  -  Dw  <*n 


fy(y  -  1)  Q)  ag 
yg'  +  y(p  “  !)“»  an 


:  12) 


yg'  +  y(p  *  Dw  an  yg'  +  y(p  -l)u  an 


(13) 


(14) 


(15) 


IV.  STABILITY  CRITERION 

The  standard  values  of  parameters  are  given  in  Table  1.  The  imaginary  part 
of  frequency  o\  and  the  wave  speed  c  are  computed  for  different  values  of  k, 
Hj/Hj,  and  H,/Hj.  Figure  3  shows  the  surface  o\  =  0  or  c  =  0  in  the  three-di¬ 
mensional  space  (k,  H2/H2,  h^/Hj).  The  surface  of  o\  =  0  is  coincident  with 
the  surface  c  =  0.  Therefore,  this  surface  divides  the  space  into  two  parts  cor- 
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TABLE  1 

The  Standard  Model  Parameters 


responding  to  a  growing  and  negative  propagation  mode,  and  a  decaying  and 
positive  propagation  mode.  Here,  the  negative  (positive)  propagation  means  that 
the  wave  propagates  in  the  direction  with  the  ice-edge  on  the  right  (left)  in  the 
Northern  Hemisphere,  and  with  the  ice-edge  on  the  left  (right)  in  the  Southern 
Hemisphere. 

Figure  3  indicates  the  following  facts,  (a)  A  negative  (positive)  propagating, 
self-excited  coastal  wave  is  growing  (decaying),  (b)  For  an  iceward  convex  ice- 
edge  (k  <  0)  the  coastally  trapped  waves  are  all  growing,  (c)  For  a  seaward  con¬ 
vex  ice-edge  (k  >  0)  the  waves  are  growing  in  a  thick-ice  region  (Hj/H2  <  0.7) 
and  decaying  in  a  thin  ice  region  (H2/H2  >  0.7). 


FIGURE  3.  Separation  of  growing  (decay¬ 
ing)  and  positively  (negatively)  propagating 
modes. 
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TABLE  2 

Some  Typical  Values  of  Growth  Rate  (Doubling  Tim  or  an  leeward  Convex  Ice-Edge 


Small  upper- 
layer  depth 
(H,/H,  =  0.01) 


Large  upper- 
layer  depth 
(H./H,  =  0.09) 


(K  = 

-0.2) 

Thin  ice 

Thick 

(H,*/H,  =  0.8) 

(H ,'/  =0.S) 

lO'V 

10  1 

(19.3  hr) 

(8  Jays) 

0.6X10‘V 

10  •' 

(134  days) 

(2  5  yr.) 

TABLE  3 

Some  Typical  Values  of  Growth  Rate  (Doubling  Time  r  Decay  Rate  (Half-life  Period*) 
for  a  Seaward  Convex  Ice-Edg  k  =0.2) 


Small  upper- 
layer  depth 
(H,/H,  =0.01) 


Large  upper- 
layer  depth 
(H./H,  =0.09) 


1 - 1 - 

Thin  ice 
(H,’/H,  =0.8) 

Thick  e 
(H =0.S) 

— 0.S6X  10~5s~' 

(1.43  days)* 

H  s'1 
(8  days) 

-0.36X10*7  s'1 
(372.2  days)* 

10  s"1 

(2  ?6  yr.) 

Very  thick  ice 

=0.1) 


0.31 6X  1 0"4  s" 
(2S.4  days) 


10'V 

(2.196  yr.) 
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FIGURES.  Isolines  of  (a)  — 0|  ( 1 0-6 s-1 ) 
for  thin  ice  (Hj/Hj  =  0.8),  (b)  logl0  ( — cr,)  (o, 
in  s'1)  for  thick  ice  (H2/Hj  =  0.5),  and  (c) 
l°9xo  < — CTi)  (0\  in  s'1)  for  very  thick  ice 
(H'2/Hj  =0.1). 
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V.  PROPERTIES  OF  THE  SELF-EXCITED  COASTAL  TRAPPED  WAVES 

Figure  4  indicates  the  isolines  of  wave  speed  with  H,/H2  and  k  at  three  dif¬ 
ferent  values  of  H2/H2  (0.8,  0.5,  0.1).  The  isoline  pattern  at  H2/H2  =  0.8  (thin 
shelf-ice)  is  quite  different  from  that  at  H2  /H2  =  0.5  and  0. 1  (thick  and  very  thick 
shelf-ice).  For  thin  shelf-ice  (Fig.  4a)  the  curvature  of  the  ice-edge  affects  the 
phase  velocity  dramatically.  It  not  only  changes  the  wave  speed,  but  also  alters 
the  direction  of  propagation.  The  coastal  waves  propagate  positively  (negatively) 
for  convex  seaward  (iceward)  shelf-ice  edge.  The  increase  in  the  absolute  value  of 
the  ice-edge  curvature  augments  the  wave  speed.  Stronger  effects  can  be  seen  in 
the  negative  area  (i.e.,  convex  iceward).  For  the  same  curvature,  the  wave  speed 
decreases  with  increasing  upper-layer  depth  H , .  However,  the  isoline  patterns  at 
H2/H2  =  0.5  and  0.1  (Fig.  4b  and  c)  are  very  similar.  The  coastal  waves  all 
propagate  negatively  regardless  of  whether  the  ice  edge  is  seaward  or  iceward 
convex.  The  wave  speed  (absolute  value)  increases  monotonically  with  decreasing 
k..  This  means  that  for  a  given  H, ,  the  coastal  wave  has  its  minimum  speed  at  k  = 
=  0.2  and  maximum  speed  at  k  =  -0.2.  For  a  given  curvature  k,  the  wave  speed 
decreases  slightly  with  the  increase  in  the  upper-layer  depth.  Figure  4c  also  shows 
that  for  very  thick  shelf-ice  (H2/H2  =0.1)  the  wave  speed  is  almost  independent 
of  H  if  H,  exceeds  some  critical  value  near  0.01  H2. 

Figure  5a  displays  the  isolines  of  the  growth  rate  —  o\  in  the  (H^Hj,  k) 
plane  for  thin  shelf-ice  (H'2/H2  =0.8).  A  zero  line  divides  the  plane  into  two 
parts:  growing  and  decaying  modes.  The  growing  (decaying)  mode  is  generated 


FIGURE  6.  Dependence  of  log10|aj|on 
H,/Hj  for  thin  shelf-ice  (Hi/H2 
■  0.8):  a)  k  ■  -0.2,  and  b)  k  -  0.2. 
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mainly  in  an  iceward  (seaward)  convex  ice-edge.  The  growth  or  decay  rate 
decreases  rapidly  with  an  increase  in  upper-layer  depth  H,.  This  is  clearly  seen 
in  Fig.  6. 

Figure  5b  and  c  shows  the  isolines  of  the  logarithmic  growth  rate  log  10  ( — o\) 
in  the  (H1/H2,  k)  plane  for  thick  (H2/H2  =  0.5)  and  very  thick  (H2/H2  =  0.1) 
shelf-ice.  The  isoline  patterns  are  very  similar  in  these  two  cases.  The  growth 
rate  ( — ct,- )  is  positive  throughout  the  whole  plane  and  depends  more  on  the  upper- 
layer  depth  than  on  the  curvature.  For  small  upper-layer  depth  H,/H2  =  0.01, 
the  growth  rate  |“-<7j)  is  around  10'6  s-1  (0.316X10-6  s_I)  in  thick  (very  thick) 
shelf-ice.  For  a  large  upper-layer  depth  H1/H2  =  0.1,  the  growth  rate  ( — C7j)  be¬ 
comes  as  small  as  10-8  s_1  in  both  the  thick  and  very  thick  shelf-ice  cases. 

Figure  6  shows  the  dependence  of  logarithmic  growth  or  decay  rate  log ,0  icr j | 
on  H,/H2  for  thin  shelf-ice  (H2/H2  =0.8).  In  the  convex  iceward  case  (k  =  -0.2), 
the  growth  rate  drops  rapidly  with  the  increase  in  the  upper-layer  depth  from 
1CT5  s"‘  (at  H/H  =  0.01)  to  0.6X1CT7  s"1  (at  Ht/H2  =  0.09).  In  the  convex 
seaward  case  (k  =  0.2),  the  damping  rate  also  decreases  rapidly  with  the  increase 
in  upper-layer  depth  from  — 0.56X10-5  s_1  (at  H,/H2  =0.01)  to  — 0.36X1  CT7  s_1 
(at  H,/H2  =0.09). 


VI.  CONCLUSIONS 

This  air-sea  coupled  model  is  intended  to  depict  the  generation  of  unstable 
coastal  water  waves  without  external  forcing  near  a  fast  shelf-ice  edge.  The  model 
shows  the  following  results: 

(a)  Negatively  (positively)  propagating  coastal  waves  are  growing  (decaying). 

(b)  For  an  iceward  convex  ice-edge,  the  coastal  waves  are  all  growing.  The 
growth  rate  decreases  rapidly  with  the  increase  in  the  upper-layer  depth  Hi  and 
the  ice  thickness.  For  a  seaward  convex  ice-edge  the  coastal  waves  are  growing  for 
thick  ice  (H2/H2  <  0.7)  and  decaying  for  thin  ice  (H2/H2  >  0.7). 

(c)  From  the  investigation  of  some  typical  values  of  growth  rate  listed  in 
Tables  2  and  3,  it  is  found  that  iceward  convex  ice-edge,  thin  ice,  and  small  upper- 
layer  depth  favor  the  generation  of  the  most  unstable  modes  of  the  coupled  air- 
sea  system,  and  seaward  convex  ice-edge,  thin  ice,  and  small  upper-layer  depth 
favor  the  damping  modes. 

(d)  The  horizontal  surface  temperature  gradient  across  the  ice-edge  also 
generates  surface  winds,  which  in  turn  produce  the  surface  water  current  and 
affect  the  vertical  displacement  of  pycnocline.  This  effect,  which  is  excluded  in 
this  paper,  should  be  taken  into  consideration  in  future  research. 
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